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ABSTRACT 

There has been growing observational evidence for the presence of condensates in the 
atmospheres and/or comet-like tails of extrasolar planets. As a result, systematic and 
homogeneous tables of dust properties are useful in order to facilitate further obser¬ 
vational and theoretical studies. In this paper we present calculations and analysis of 
non-isotropic phase functions, asymmetry parameter (mean cosine of the scattering 
angle), absorption and scattering opacities, single scattering albedos, equilibrium tem¬ 
peratures, and radiative accelerations of dust grains relevant for extrasolar planets. 
Our assumptions include spherical grain shape, Deirmendjian particle size distribution, 
and Mie theory. We consider several species: corundum/alumina, perovskite, olivines 
with 0% and 50% iron content, pyroxenes with 0%, 20% and 60% iron content, pure 
iron, carbon at two different temperatures, water ice, liquid water, and ammonia. The 
presented tables cover the wavelength range of 0.2 to 500 micron and modal particle 
radii from 0.01 micron to 100 micron. Equilibrium temperatures and radiative acceler¬ 
ations assume irradiation by a non-black-body source of light with temperatures from 
7000K to TOOK seen at solid angles from 27r to 10“^ sr. The tables are provided to the 
community together with a simple code which allows for an optional, hnite, angular 
dimension of the source of light (star) in the phase function. 

Key words: astronomical data bases: miscellaneous - opacity - scattering - planets 
and satellites: atmospheres - circumstehar matter. 


1 INTRODUCTION 

Many classes of astronomical objects have low enough tem¬ 
perature and high enough density that grains of conden¬ 
sates can be formed. Such grains are usually called “dust”, 
although some authors use the more generic term “conden¬ 
sates”. These objects include interstellar or interplanetary 
medium, molecular clouds, the outer regions of active galac¬ 
tic nuclei discs, outer atmospheres of cool stars (e.g., Mira- 
type supergiants), supernova remnants, brown dwarfs, ex- 
trasolar planets, circumstehar, protoplanet ary, and debris 
discs, comets, and others. 

The optical properties of condensates may not only in¬ 
fluence but fully govern the emerging spectrum and struc¬ 
ture of the object. Dust can absorb the impinging radia¬ 
tion and convert it directly into heating of the grains. This 
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process is called ’absorption’ or ’true absorption’ to em¬ 
phasize that the photon is thermalized. It is quantihed by 
the absorption opacity. Dust can also scatter radiation in a 
process called ’scattering’ without being heated. This pro¬ 
cess is characterized by the scattering opacity. Furthermore, 
scattering can be highly asymmetric, a property that is de¬ 
scribed by means of the phase function, which depends on 
the scattering angle (the deflection angle from the original 
direction of the impinging radiation). The most prominent 
feature is a strong forward scattering for large values of 
X = I'KajX where a is particle size (radius) and A is wave¬ 
length. Finally, condensation affects the chemical composi¬ 
tion of the object. It removes the condensed elements from 
the gas phase within the dust cloud and, in the atmosphere, 
also from the region above the clouds due to the rain-out. 

This paper focuses mainly on extrasolar planets. Ex¬ 
trasolar planets (exoplanets) are found at a wide range of 
orbital periods and/or distances from their host stars, rang- 
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ing from a fraction of a day to many years. Their hosts have 
spectral types from A-stars to M-stars and down to brown 
dwarfs. As a result, depending on the amount of stellar ir¬ 
radiation, the atmospheres of exoplanets may have temper¬ 
atures from a few thousand K to about one hundred K. At 
such temperatures condensates are expected to be present 
in their atmospheres or circum-planetary material. Conden¬ 
sates are usually confined to clouds at a certain depth or re¬ 
gion according to their condensation temperatures (see e.g. 
Ackerman fc Marley||20QT] [Allard et al.| 2001| |Tsuji|[2002 


which moves the clouds and the L to T spectral transition 
to lower effective temperatures than in brown dwarfs. 

Forward scattering on dust can be particularly impor¬ 
tant in eclipsing systems during the primary eclipse, when 
a cool dusty object is in front of the main source of light 
( Budaj [2011 ). Recently a few such systems with dusty disks 
were discovered ( Dong et al.||2014 Meng et al. 2014 Rat- 
tenbury et al.|[2015 ). In extrasolar planets, this situation is 


Lodders|2003| j Sanchez-Lavega, Perez-Hoyos &: Huesoj2004 


exactly encountered near the transit, de Kok & St am (2012) 


Burrows, Sudarsky Hubeny12006||Lodders k, Fegley|2006 

Helling et al.||2008||Morley et al.||2012||2014| ). In highly ir- 


studied how forward scattering affects the transmission spec¬ 
troscopy of exoplanet ary atmospheres. However, dust clouds 
associated with planets are not necessarily confined to the 


radiated atmospheres, or at large depths of cooler exoplan¬ 
ets, only the most refractory dust species (Ca-Al-Ti bearing 
compounds) can survive. They are overlaid by the optically 
opaque clouds of iron and various silicates. Higher up in 
the atmosphere, or at cooler temperatures, finer sulfide and 
alkali halide clouds may be present. Water clouds become 
dominant at lower temperatures and, at even lower tempera¬ 
tures, ammonia clouds are expected. An entire classification 


atmosphere of the planet. Rappaport et al. (2012) discovered 


a transiting disintegrating exoplanet KIC 12557548b with a 
comet-like dusty tail extending well beyond its Hill radius 
and featuring a pre-transit brightening. Recent studies have 
shown that pre-transit brightening is due to strong forward 
scattering and can be used, together with the color depen¬ 
dence of the transit depth, to constrain the particle size of 
dust grains in the comet-like tail, which was found to be of 


on their spectra. 

The above-mentioned dust phase function should not 
be confused with the planet phase function. The later de¬ 
scribes the light curve of a planet as a function of orbital 
phase and is modulated by the dust phase function. There 
are numerous theoretical studies and numerical models of 
the light curves, spectra, planet phase functions, and albedos 


scheme for exoplanets was proposed by [Sudarsky, Burrows[ 

ski et al.[[2015[) (see also Brogi et al.[[2012 van Werkhoven 

& Pinto (2000) based on the expected effects of such clouds 

et al.||2014|). Two other objects of this kind (KOI-2700b, 


EPIC201637175B) have been discovered already by Rappa- 
jport et ^ ( |2014| and Sanchis-Ojeda et al. ( 2015| ). The ob¬ 
served tail lengths in these objects are consistent with corun¬ 


dum or iron-rich silicate dust grains (van Lieshout, Min & 
Dominil^|2014 ). Finally, scattering on dust in the Earth at¬ 


mosphere may have important effects on lunar eclipses or on 
the Earth transmission spectr um (jKocifaj |1994| Kocifaj & 


in extrasolar planets including the effect of clouds (Marley 
et al.|199'^|Seager, Whitney 5asselov|2Q00[|Sudarsky et al. 


2005 Cahoy, Marley Eortney|2010| Kane L, Gelino|2011| ). 

A convenient analytical model for albedos and phase curves 
taking into account non-isotropic scattering was presented 


Horvath|2005||Palle et ah' 2011[|Garcfa Muhoz Palle|2011| 

Garcfa Munoz et al.|2011 ), as well as on the light curve and 

spectra of Earth ( Tinetti et al.||2006| Stam||2008 ) . A recent 
summary of the cloud and haze formation in exoplanet at- 


Madhusudhan & Burrows (2012) and Heng, Mendonga & 


Lee (2014). It turned out that hot-Jupiters have very small 


albedos in comparison with the Solar System planets and it 
is quite difficult to detect them in reflected light in the opti¬ 
cal region due to strong absorption in the Na, K resonance 
lines and/or absence of highly reflective condensate grains 
( Sudarsky, Burrows Pinto|2000 Rowe et al.|2008 Cowan 


Agol|2011| . Nevertheless, recent observations support the 
presence of clouds in their atmospheres. Pont et al. (2013) 
and Lee et al. (2014) state that transmission spectra of 


mospheres can be found in Marley et al. (2013). 

Radiative transfer calculations in planetary atmo¬ 
spheres are very CPU-time demanding. As a result, in many 
applications researchers use simplifying assumptions such 
as isotropic scattering or analytical Henyey-Creenstein dust 
phase functions given that more exact calculations require 
application of Mie theory or the use of additional Mie- 
scattering codes. In such cases, precomputed dust opacities 
and phase functions may be very convenient and efficient 


HD 189733 are consistent with the presence of condensates in 


the atmosphere of this planet. Sing et al. (2013) argues that 


aerosols, especially corundum, may explain the transmis¬ 
sion spectra of another hot-Jupiter, WASP-12b. |Wakeford fcj 
[Sing ( |2Q15| examined the effect of various grain species and 
their particle sizes in the atmosphere of a typical hot Jupiter 
(HD 189733b) on the transmission spectrum. They predicted 
strong enstatite and perovskite features in the infrared re¬ 
gion. Also, the atmosphere of super Earth CJ 1214b must 
contain clouds to be consistent with the observed transmis¬ 
sion data ( Kreidberg et al.||2014 Morley et al.||2013| [Howe 


( Hubeny Lanz|1995 Burrows, Sudarsky Hubeny|2Q0'6 ). 
The aim of this paper is to facilitate further research and 
provide the community with homogeneous and systematic 
tables of dust phase functions, opacities, single scattering 
albedos, grain equilibrium temperatures, and radiative ac¬ 
celerations relevant mainly for the field of extrasolar planets 
or circum-planetary (circum-stellar) material. These tables 
can be incorporated into more complex cloud models, or 
find additional applications in brown dwarfs, protoplanet ary 
discs, comets, and interacting binary stars. 

The supplied tables may also be beneficial for anal¬ 
ysis and interpretation of observations from future space 
missions. Large area planet transit surveys such as TESS 


fc Burrows||2012 ). Madhusudhan, Burrows fc CurrieJ (2011) 

concluded that directly-imaged HR 8799 planets have thick 
clouds (thicker than those required to explain data for typi¬ 
cal L and T dwarfs). On the other hand, Marley et al. (2012) 
suggested that it is because exoplanets have lower gravities 


(Ricker et al. 2014) and PLATO (Rauer et al. 2014) will 


discover many exoplanets that will be followed-up by other 
large telescopes and their atmospheres, albedos, and phase 


curves studied (Deming et al. 2009). Superb photometric 


precision of CHEOPS (Broeg et al 


2013) will detect exo¬ 


planetary atmospheres and study day-night heat transfer in 
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hot-Jupiters, and their albedos. JWST ( Gardner et al.|2006 
Beichman et al. 2014) will perform transit or occultation 


spectroscopy and direct imaging of planets which will enable 
chemical composition analysis of the planetary atmospheres. 
It will study dust in all sorts of environments including our 
own solar system or proto-planetary disks. Another mission, 
ECHO, dedicated to study the composition of atmospheres 


of transiting exoplanets was proposed (Tinetti et al. 2012 


2015). 


Since there are currently no observational constraints on 
the shape, porosity, or aggregate structure of grains in these 
environments, we assume homogeneous spherical particles 


and Mie theory. Recently, Cuzzi, Estrada & Davis (2014) 


proposed an opacity model which takes into account porous 


aggregates, Mordasini (2014) developed an analytical model 
for grain opacity in the atmospheres of forming planets, and 


Ormel (2014) presented a method to include the evolution of 


grain size and opacity into the structure of protoplanetary 
atmospheres. The above-mentioned features are not contem¬ 
plated in the present paper and it was not our intention to 
calculate cloud or atmosphere models for exoplanets. 

This paper is organized as follows. Section describes 
the source of the refractive index. Section [^describes the pa¬ 
rameter space covered by the calculations. Sections de¬ 
scribe the calculations and properties of the phase functions, 
opacities, albedos, equilibrium temperatures, and radiative 
accelerations of each species. Section provides a detailed 
description of the tables, their format and units. Sections 
contains our summary and conclusions. 


2 COMPLEX REFRACTIVE INDEX 

The optical properties (cross-sections for scattering, ab¬ 
sorption, and phase functions) of spherical dust particles 
composed of a homogeneous material with a wavelength- 
dependent refractive index can be calculated using the con¬ 
ventional Mie theory. Here we describe the sources of the 
refractive index. The Heidelberg - Jena - St.Petersburg - 
Database of Optical Constants is a very convenient starting 
point ( Henning et al.|19^ Jager et ^|2Q03b ). 

Corundum. It is the most common form of crystalline 
alumina (AI 2 O 3 ), which is a alumina. It is one of the most 
refractory condensates that is expected to condense from 
a gas of solar chemical composition. Corundum is an op¬ 
tically anisotropic (birefringent) material. We adopted re¬ 
cent IR-mm measurements of its refractive index by |Zeidler^ 
Posch &: Mutschk^ ( |2013| ). These are in situ, high temper¬ 
ature measurements, in the range of 300-928 K for both 
ordinary and extraordinary rays. The optical properties of 
corundum were found to be temperature dependent. Corun¬ 
dum is most important for exoplanets at high temperatures 
since at lower temperatures its effects are superseded by 
more opaque and less refractory species. For this reason, we 
adopted the measurements for the highest temperature con¬ 
sidered in the above mentioned study (928 K). We calculated 
the optical properties (cross-sections, opacities, phase func¬ 
tions, asymmetry parameter) independently for both ordi¬ 
nary and extraordinary rays. Then we took a weighted sum 
of them with weights of 1/3 and 2/3 for extraordinary and 
ordinary rays, respectively, as discussed by |Zeidler, Posch 
Mutschke (2013) to account for an averaged crystal orienta¬ 


tion. There are several types of alumina which may differ in 
optical properties. Therefore we considered also 7 alumina 
and its complex refractive index was taken from |Koike et ah'] 
(1995). Their measurements are at room temperature and 


cover a broad range of wavelengths from UV, optical to far 
IR. If not otherwise specified, in the next, we will refer to 7 
alumina as alumina and a alumina as corundum. 

Perovskite. Perovskite is calcium titanium oxide (cal¬ 
cium titanate) with the chemical formula CaTiOa. Similar 
to corundum, it is one of the most refractory species and the 
highest temperature condensates from a gas of solar chemi¬ 
cal composition. It may form clouds in exoplanets and brown 
dwarfs ( Fodders Fegley|2006 ) which might be detected in 
the transmission spectrum in the infrared region with instru¬ 


ments like JWST (Wakeford & Sing 2015). Formation of per¬ 
ovskite is interesting also because it removes titanium and 
thus also TiO from the gas phase. TiO gas is a major opacity 
source in the optical region of M dwarfs. If present, it would 
be a strong absorber also in exoplanets, causing temperature 
inversions - stratospheres (Hubeny, Burrows & Sudarsky 


2003 Knutson et al. [20(38 Burrows, Budaj &: Hubeny|2008 ^ 
with huge impact on the atmospheric structure and spectra. 
The complex refractive index of perovskite was taken from 


Posch et al. (2003). They made reflectance spectroscopy of 


a natural perovskite crystal and used Lorentz oscillator fits 
to determine the optimum parameters, such as frequencies, 
strengths, or damping constants. Lorentz oscillator param¬ 
eters were used to compute optical constants n and /c, that 
in turn predetermine small particle spectra. Unfortunately, 
the dielectric constants obtained this way do not cover the 
whole spectral region and span only wavelengths longer than 
2 microns. 

Pyroxenes. The complex refractive index of pyroxenes 
was taken from Dorschner et al. (1995), who determined 


the optical constants from reflectance and transmittance 
measurements. While the reflectance data were available 
for wavelengths ranging from 0.2 jjim to 500 fim with 
Kramers-Kronig analysis in 8-80 /xm, the transmittance 
analysis was possible only at wavelengths below 8 (im and 
then above 100 microns. The quality of data sets was guar¬ 
anteed by homogeneity tests made for the sample dur¬ 
ing lab experiments. The optical properties of pyroxene 
(MgxFei-xSiOs) in the optical region are quite sensitive to 
the amount of iron 1 — x in the mineral. Therefore, we con¬ 
sidered three kinds of pyroxenes: enstatite (iron free, x = 1 ), 
pyroxene with 20 % of iron (x = 0 . 8 ), and pyroxene with 60% 
of Mg atoms replaced by iron (x = 0.4). 

Olivines. Forsterite (Mg 2 Si 04 ) was adopted as an iron- 
free representative of the olivine family of magnesium sili¬ 


cates. Its optical data were taken from Jager et al. (2003a), 


who used sol-gel technique to produce amorphous magne¬ 
sium silicates and followed with homogeneity and purity 
tests based on TEM analysis combined with EDX. We con¬ 
sidered also an iron enriched olivine (Mg 2 ^Fe 2 - 22 /Si 04 ) with 
half-half of Mg-Fe atoms {y = 0.5) and the data were taken 
from Dorschner et al. 

Iron. In contrast to silicates, iron-like particles show dif¬ 
ferent optical characteristics as discussed by | Pollack et al.| 
(1994), where the real and imaginary values of complex re¬ 


fractive index were published for wavelengths ranging from 
0.1 to 10^ iim. However, for pure iron we have used two data 


sources: Johnson & Christy (1974) for A < 0.7 fim and Ordal 
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et al. (1988) for A > 0.7 fim. In both cases, the experimental 


data were collected at room temperature. For the metals, the 
Drude model - the classical phenomenological approach for 
the conductivity - is traditionally used to explain the trans¬ 
port properties of electrons. However, the Drude model is 
not well suited for the near and mid infrared region, thus it 
may potentially be inaccurate when applied to A > 0.7 jim. 
Nevertheless, Ordal et al. (1988) concluded that, within the 
experimental error, the non-resonant cavity data can be fit 
with plasma frequency 32,900 cm“^ and damping frequency 
213 cm“^. Combining their measurements with Kramers- 
Kronig analysis, the frequency-dependent refractive indices 
were computed and provided along with the reflectance data. 

Carbon. Carbon, mainly in the form of graphite parti¬ 
cles, is an important ingredient of interstellar dust. It has 
not been detected in extrasolar planets, yet. Nevertheless, 
we included carbon dust particles into our tables since car¬ 
bon grains are potentially abundant refractory species and 
might be present in circumstellar material, carbon-rich en¬ 
vironments, or disintegrating exoplanets. The complex re¬ 
fractive index of carbon dust may be quite sensitive to the 
temperature due to the presence of free charge carriers at 
higher temperatures as shown by |Jager, Mutschke Hen-| 


nmg 


(1998). These authors measured the optical properties 
of amorphous carbon for a set of temperatures ranging from 
400 C to 1000 C for a broad wavelength range from 0.2 to 
about 500 microns. As the temperature increased, an in¬ 
creased carbonization and graphitization of the sample was 
observed. We adopted their measurements at the highest 
and lowest temperatures. 

Water. A new update of the water ice optical constants 
was published recently by Warren & Brandt (2008), who re¬ 
vised the imaginary part of the complex refractive index in 
the visible spectrum. The values they tabulated are appro¬ 
priate for the hexagonal crystal form of ordinary ice or frozen 
water at temperatures not far below the melting point (the 
nominal temperature was 266 K). Regarding liquid water 
drops, the optical data published by the International As¬ 
sociation for the Properties of Water and Steam (lAPWS) 
are probably the most accurate, but are available only for 
the real part of the refractive index and valid only for a 
limited spectral range. Therefore we have used Segelstein 
data ( Segelsteln||1981 ), which cover a much wider range of 
wavelengths for both real and imaginary values of the refrac¬ 
tive index. Segelstein data are well accepted by the optical 
community and thus often used in simulating the optical 
properties of liquid water in various environments. 

Ammonia. For ammonia (NH 3 ) at temperature ranging 
from 77 K to 88 K we used [Martonchik, Orton k, Appleb;^ 
(1984), which is the more comprehensive review of the op¬ 
tical constants of NH 3 ice. However, we stress that optical 
constants for solid ammonia in the cubic phase may differ 
significantly from those in other phases - as shown by |Fer-| 


raro. Sill & Fink (1980). 


3 PARAMETER SPACE - GRID. 

Particle size. We assumed homogeneous spherical particles 
with a poly-dispersed Deirmendjian particle size (radius) 
distribution ( Deirmendjian|19^ ), Ur, given by: 



Particle radius [mic] 


Figure 1. Deirmendjian particle size distribution with a modal 
particle size of 1 micron. 



where rc is the critical or modal particle radius, where the 
function has its maximum. We truncate the size distribution 
at 0.22 rc for small particles and at 6rc for large particles. 
Particular attention is paid to the extended tail of larger par¬ 
ticles, since their contribution may be significant e.g. to the 
Rayleigh scattering or to the forward scattering peak. This 
particle size distribution is shown in Fig[^ We define the 
distribution on a very fine integration step in order to aver¬ 
age the phase functions and to suppress the ripple structure 
that would appear in the phase function of spherical mono- 
dispersed particles. For large particle sizes it is necessary to 
take into account as many as 1600 points in the particle size 
distribution. We calculate the dust optical properties for 21 
particle size distributions with different modal radii ranging 
from 0.01 to 100 micron with an equidistant step of 0.2 dex 
in the logarithm of rc. Such particle size distributions can 
be expected for example in planetary atmospheres. Note, 
however, that although the range of our particle size distri¬ 
bution (see Fig[^ spans more than 1 dex, the typical full 
width at half maximum of such a distribution is relatively 
narrow, about 0.5 dex. Consequently, our results obtained 
assuming a distribution with a certain modal radius might 
be used also as a representative values for particles with radii 
equal to the modal particle radius and might be assembled 
to construct optional, broader, particle size distributions if 
necessary. 

Wavelength coverage. Optical properties are calculated 
for 400 frequency points covering the region 0.2 to 500 mi¬ 
cron, and spaced equidistantly in logarithm of wavelength. 
This wavelength range covers most of the spectral energy 
distribution (SED) of stars with effective temperatures lower 
than 8000 K. Not always the entire wavelength region is 
covered by the measurements of the refraction index (see the 
Tablefor summary). 

Scattering angles. Phase functions are calculated for 65 
angles ranging from 0 to 180 degrees. Since the scattering 
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Table 1. Summary of the refraction index measurements adopted 
in this work. 


species 

A 

N 

pg 

reference 

alumina 

0.2-400 

617 

2.9 

K 

corundum 

6.7-le4 

1500 

2.9 

Z 

perovskite 

2.0-5843 

2000 

4.1 

P 

iron 

0.19-1.9 

49 


J&C 

iron 

0.62-286 

53 

7.874 

O 

forsterite 

0.20-949 

4847 

3.0 

JOS 

olivine50 

0.20-500 

109 

3.0 

D 

enstatite 

0.20-500 

109 

3.0 

D 

pyroxene20 

0.20-500 

109 

3.0 

D 

pyroxene60 

0.20-500 

109 

3.0 

D 

carbonlOOO 

0.20-794 

73 

1.988 

J98 

carbon0400 

0.20-518 

69 

1.435 

J98 

water ice 

0.04-2e6 

486 

1.0 

W&B 

water liq 

0.01-le7 

1261 

1.0 

S 

ammonia 

0.14-200 

298 

0.88 

M 


Type of grains, A -wavelength range [/rm], N -number of wave¬ 
length measurements, p® -grain density [gcm~^], reference: K - 


Koike et al. 

(|1995|), Z - |Zeidler, Posch & Mutschke 

(|2013|), P - 

Posch et al. 

(|2003|), J&C -|Johnson & Christy] ( 

1974 

1 , 0 -|Ordal 

et al. 

(|198g 

I, JOS - |Jager et al.| (|2003a 

, D - 

Dorschner et al. 

(|1995 

1 , J98 -|Jager, Mutschke & Henning 

1 ( 19 ^ 

, W&B -1 Warren 

& Brandt 1 ( 

2008|), S - |Segelstein| (|1981|), M - |Martonchik, Orton 

& Appleby 

|T9^. 


is not isotropic, and may have a very strong and narrow 
forward scattering peak, we adopt a much finer step in the 
forward and backward scattering regimes. Angles are sym¬ 
metric with respect to 90 degrees. 

Dust species. Optical properties of different condensates 
are calculated on the same grid of 21 particle size distribu¬ 
tions and 400 wavelengths. They are: corundum/alumina, 
perovskite, pure iron, forsterite as an iron-free representa¬ 
tive of the olivine family of silicates, olivine with 50% of 
iron, enstatite - an iron-free representative of the pyroxene 
family of silicates, pyroxenes with 20% and 60% of iron, car¬ 
bon at two different temperatures 1000 C and 400C, water 
ice, water liquid, and ammonia. 

Star temperatures. In addition to inherent dust param¬ 
eters (dust chemical composition, particle size) some grain 
properties (equilibrium grain temperature, radiative accel¬ 
eration) depend also on the properties of the irradiating ob¬ 
ject. The most important property of the light source is its 
SED and its angular dimension (i.e. its size and distance). 
The SED on the surface of the source is mainly a func¬ 
tion of its effective temperature and, to a lesser extent, also 
of the surface gravity, chemical composition, microturbu¬ 
lence and/or other free parameters. In case of theoretical 
models the output SED depends also on various simplyfy- 
ing assumptions and/or input data. Eor this reason, we cal¬ 
culated equilibrium grain temperatures for 13 star/brown 
dwarf effective temperatures: 7000, 5800, 5000, 4500, 4000, 
3500, 3000, 2500, 2000, 1600, 1200, 900, and 700 K. Eor 
objects with temperatures ^ 2500 K we assumed the sur¬ 
face gravity \ogg = 5.0 (cgs), while for hotter objects we 
assumed log^ = 4.5 (cgs). The source is further assumed 
to have a solar chemical composition ( Asplund et al.|[2009 


Caffau et al.|2011 ). We use a single homogeneous grid of syn¬ 
thetic spectra based on atmospheric model calculations, BT- 


Settl (CIEIST2011 -2015, CIEIST2011) ( |Ahard et alT1[2Q03 
Baraffe et al.|201^ , which covers the entire parameter space. 
These models are ID, static, with spherical symmetry, in lo¬ 
cal thermodynamical equilibrium, with convection and take 
into account atomic, molecular and dust opacities. The spec¬ 
tra list energy hux, Fx, on the surface of the object in 
erg s~^cm~^cm~^ from 0.001 to 1000 micron. This can be 
easily converted to, in units of erg s~^cm~^Hz~^ via 
Fi, — FxX^/c. We recall that the effective temperature, Teff, 
of the object is, by dehnition, the temperature of the black- 
body with the same total energy output i.e.: 


/ 


Fjydiy = crTeff = TT / Bjy{Tefi)diy, 


/■ 


( 2 ) 


where a is the Stefan-Boltzmann constant and Bjy is the 
Planck function. 

Solid angles. The dust temperature depends also on the 
solid angle subtended by the source of light. We considered 
15 values, roughly logarithmically spaced, ranging from 27r 
up to 10sr. Radiative accelerations are calculated for the 
same stellar effective temperatures and particle sizes. Their 
dependence on the solid angle becomes insignihcant at a 
relatively small distance from the source and thus was not 
included. 


4 CALCULATIONS 
4.1 Opacities 

We use a widely available Mie scattering code CALLBH- 
MIE adapted by B.T Draine. It calls iteratively the Bohren- 
Huffman Mie scattering subroutine BHMIE, |Bohren Huff-| 
(1983), Appendix A. It provides Qa, Qs, and Qe, which 


are efficiency factors for absorption, scattering, and extinc¬ 
tion, respectively, for particles with radius r. They are re¬ 
lated to the Ca and Cs, the absorption and scattering cross- 
sections, via 


Ca = QaTrr^, Cs = QsTVr^, Qe = Qa + Qs 


(3) 


It is assumed that the dust grains are composed of homo¬ 
geneous material, have spherical shape, and are placed in 
an environment with the complex refractive index of 1-Oi. 
The code is written in EORTRAN 77. It was modihed to 
calculate the average cross-sections for absorption Ca and 
scattering Cs for an ensemble of particles with an optional 
size distribution Ur- 


Ca=- 

n 




CaTirdr 


- 

n 


d 


CsUrdr, 


where 

n = J Urdr. 


(4) 


(5) 


Subsequently, the absorption and scattering opacities 


of condensates in units of cm ^ are given by 




_dust /-Y _ 

cr,. = GsU. 


( 6 ) 


Eor practical applications it is more convenient to use the 
opacities per gram of dust material, given by 


dust dust / 

^i/,p ^1/ / Pd 5 


dust dust / 

/ Pd 


( 7 ) 
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where pd is the density of dust made of particular species 
and 1 are absorption and scattering opacities, re¬ 
spectively, per gram of dust material in units of cw? g~^. 
The dust density can be expressed as 


pd = 


J MrUrdr = p^ J V^Tirdr — p^V^n — M^n 


(8) 


where Mr is mass of a dust grain of radius r. are the 

mean volume and mass of a dust grain, respectively, and 

Ms = — J V^Tirdr = J r^Tirdr. (9) 

In the above expression, p^ is the (constant) density of a 
grain of dust. Substituting equations § and into equa¬ 
tion 0 one can obtain the following relation between opac¬ 
ities and cross-sections: 


dust _ 


Ms 


( 10 ) 


( 11 ) 


dust Cs 

MS 

Using these opacities, the monochromatic optical depth 
along the line of sight z is then given by: 


/ 


= - / p4z)[K^^\z) + atTid]dz. 


( 12 ) 


In the following, the sum of the absorption and scattering 
opacities will be referred to as a total opacity. Table also 
lists the densities of dust grains, p®, which we assumed in 
the opacity calculations. Note that the density of corun¬ 
dum is Ag cm~^ but since corundum is most interesting 
at high temperatures when other species are still in the gas 
phase, we adopted the density of corundum close to its melt¬ 
ing point which is significantly lower, about 2.9 g cm~^. On 
the contrary, the density of ammonia ice rises significantly 
with temperature, but reaches a plateau be tween 60-100K 


at about 0.88 gem (Satorre et al. 2013), which is the 


value we adopted. The density of amorphous carbon was 
also found to increase significantly with temperature due 
to increasing graphitization. We adopted the densities from 


Jager, Mutschke & Henning (1998) for the temperatures of 


400 C and 1000 C which correspond to those adopted for 
the complex refractive index. 


4.2 Phase functions 

We used the same code to calculate the Sn element of the 
scattering matrix for uniformly sized particles with radius 
r. Sii describes the angular distribution of the scattered 
light assuming that the incident light is unpolarized. It is 
convenient to normalize this angular distribution such that 
its integral over the whole solid angle du — sud[a)dad(t) is dvr. 
Such distribution is referred to as a phase function, p((a,r), 
where a is the scattering angle which measures the deflection 
from the original direction of the photon: 

/ n27T niT 

p{a^r)du — / / p{a^r)sui[a)dad(j) — A'k. (13) 

Jq Jo 

The phase function of the population of particles, p{a), with 
the size distribution Ur, can be calculated in the following 
way: 


p{a) = J p{a,r)Csnrdr/{Csn) (14) 

Since CsTir does not depend on direction, the averaged phase 
function p has the same normalization. 

The mean cosine of the scattering angle g, also known as 
the asymmetry parameter, is calculated from the normalized 
poly-dispersed phase function as: 


9 = 


J p{a) cos{a)du / 


/ 


p{a)duj. 


(15) 


Once the opacities are known, the thermal emissivity of 
condensates (energy per unit time, frequency, volume, and 
solid angle) associated with its true absorption can be cal¬ 
culated as: 


th.dust 

Cl. 


Bj,(T)pdK 


dust 


(16) 


where Bjp is the Planck function. The angle-dependent scat¬ 
tering emissivity of condensates into a direction cto can be 
calculated from the scattering opacity and phase function 
via the following expression: 

^sc,dust(^o) = JJ p[a{ao,ai),r]CsnrL[0{ai,(l))]drdu/{47T){17) 
or 

^sc,dust(^o) = J p(ct)pda^,p®^/i.dcd/(47r), (18) 

where ljp[0{ai, cj))] is the specific intensity of radiation com¬ 
ing from the direction (cti,^) with solid angle duj = 
sm{ai)daid(l), see Fig. It is usually determined by solv¬ 
ing the transfer equation as a function of angle 0 = 0(ai, 0) 
with respect to the normal to the surface of the object. In 
the following we will write the specific intensity simply as 
Ijp. One can define another averaged phase function. Pda 
such that: 


n / N fpia)eduj j p{a)I^duj 

= jKd^ = 47rJ. ’ 

where, is the mean intensity given by: 

Ju = j Pdcd/(47r). 


(19) 


( 20 ) 


This phase function is an average over the source of light. 
Since intensity depends on the angle, equation ( p^ does not 
guarantee that Pda is normalized, contrary to e.g. p{a^r). 
However, using Pda, the scattering emission can be ex¬ 
pressed in a very simple way as: 

^sc,dust(^o) = PDA(ao)pd^T^,p^Vi.. (21) 


If we consider a distinct source of light, such as a star for 
example, and a medium surrounding the source and our dust 
grain, which is optically thin, or if we are assuming single 
scattering events, then Pda is an average of the phase func¬ 
tion over the surface of the star and ao can be conveniently 
measured from the ray originating at the centre of the star. 
If the star is far or if it is small, then Pda equals p. How¬ 
ever, if the star has a non negligible angular dimension on 
the sky, compared to the characteristic changes in the phase 
function, then one has to take its angular dimension into ac¬ 
count. The phase function of large grains or at short wave¬ 
lengths usually has a very strong forward-scattering peak, 
which can be much sharper than the stellar disc as seen e.g. 
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from a typical hot-Jupiter or a close in extrasolar planet. To 
take this effect into account one has to split the stellar disc 
into elementary surfaces and integrate the phase function 


over the disc as in the equation (19). In doing so we assume 


a quadratic limb darkening of the stellar surface: 
lu = Ji/(0)[1 — ui{l — cosO) — U 2{1 — cos^)^], 


( 22 ) 


where li^(0) is intensity perpendicular to the surface of the 
source and 0 is angle between the line of sight and a normal 
to the surface. 

One has to be cautious and calculate the phase func¬ 
tions with a very fine step near the zero angle because of 
the strong forward scattering and, consequently, the disc 
averaged phase function with a very fine step near the edge 
of the stellar disc. An illustration of such phase functions, 
which take into account the hnite dimension of the source 
of light, is depicted in Fig. for WASPlOSb. For this ex¬ 
ample we assumed a stellar radius of W = 1A4Rq and a 
planet-to-star separation of r 4.27Rq ( Gillon et al.||2014 


Southworth et a'r]|2015 ). The quadratic limb darkening co¬ 


efficients ui = 0.36,1^2 = 0.31 were taken from [Magic et aL\ 
(2015). This results in an angular radius of the stellar disc of 


WASP103 of 19.7 degrees i.e. the whole star would appear 
on the planet’s sky spanning almost 40 degrees in diameter! 
One can see that the finite dimension of the source of fight 
is important for larger particles or shorter wavelengths, but 
has almost no effect for smaller particles and/or longer wave¬ 
length. The Fortran90 code which reads our phase function 
tables and calculates such disc averaged phase functions Pda 
is provided together with the tables. As mentioned above, 
this code and precalculated Pda phase functions are mainly 
useful in the case that a dust grain ’sees’ a distinct source of 
fight and/or in the optically thin dust regime. However, they 
might be utilized also within the planetary atmospheres, in 
which case, one could use Pda for the first scattering event 
and p for subsequent scattering events in, e.g., Monte Carlo 
radiative transfer simulations. 


4.3 Albedos 

We calculate and tabulate both the absorption and scat¬ 
tering opacity as a function of wavelength for each modal 
particle radius. These tables and quantities can be used eas¬ 
ily to calculate the single scattering albedos, vu, for each 
wavelength and modal particle radius via the expression: 

^dust 

^dust —dust ■ 

hjiy U u 



Figure 3. Comparison of phase functions assuming a point 
source of light (solid lines) versus the phase functions that take 
into account the finite dimension of the stellar disc (dotted lines). 
Example is for enstatite at 600 nm for different dust particle radii. 
Dotted line for 0.1 micron particles almost coincides with the solid 
line. The vertical line illustrates the angular radius of the stellar 
disc of WASP 103 as seen from the planet WASPlOSb ( [Gillon et al.| 
[2014[ [^uthworth et al.|2015[ >. 


If not mentioned otherwise, in the text below, by albedo we 
mean the single scattering albedo. 


4.4 Grain equilibrium temperatures 


Electromagnetic interaction with a particle results in heat¬ 
ing that is ruled by the imaginary part of the complex re¬ 
fractive index rrn of the particle. The larger mi, the more 
fight (or electromagnetic radiation) is absorbed and conse¬ 
quently transformed to other forms of energy - most typi¬ 
cally heat. Let us assume that a particle with an absorption 
cross-section Ca is irradiated by a radiation field with an 
intensity R. The energy absorbed by the particle per unit 
time is: 


dE+ 

dt 


CaRdudu — 477 


/ 


CaJudiy^ 


(24) 


where cc is solid angle. Only angles subtended by the source 
(star) count since R is zero elsewhere. Ju is the mean inten¬ 
sity given by equation ( [^ . We will assume that this energy 
is balanced by the thermal radiation of the particle at a rate: 


CaB.{T)dv, (25) 

where Bjy(T) is the Planck function and T is the tempera¬ 
ture of the grain. Multiplying both equations by the particle 
size distribution, integrating over the particle size, and di¬ 
viding by 477 one recovers the radiative equilibrium equation 
in terms of opacities: 


= 4,/ 


/ 


= j 4'^^^B4T)diy. 


(26) 
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It is common practice to assume a black body radiation for 
the star and neglect its limb darkening, in which case the 
mean intensity can be expressed as: 


J. = 


(27) 


In the above expression, T* is the effective temperature of 
the star and uj is the solid angle subtended by the star: 


(jO — 2tt \ 1 — 



(28) 


where Ri, is radius of the star and d is the distance to the 
star centre. However, radiation emerging from stars or brown 
dwarfs may have significant departures from the black-body 
radiation, especially in cooler objects with strong molecular 
and dust opacities. Brown dwarf spectra are dominated by 
the IR water bands which affect their spectral energy dis¬ 
tribution, redistribute the flux towards shorter wavelengths 
and interfere with our dust absorption. That is why, in the 
calculations of the left-hand side of equation (26), we use 


the fluxes, iT, from the model atmospheres rather than the 
Planck function. In this case, the following equation pro¬ 
vides a much better estimate of the mean intensity (recall 
equation §): 


J. = 


F^{T*) 

47r TT 


(29) 


Equation (26) can be solved iteratively to obtain the 


grain temperature as a function of stellar temperature and 
solid angle. Numerical integration of the equations was per¬ 
formed using a simple trapezoidal rule. Particles of different 
sizes will have different temperatures as a result of their dif¬ 
ferent absorption opacities. However, since our particle size 
distribution is relatively narrow, we assume that particles 
within the same distribution have the same temperature. 
Consequently, we calculate and tabulate grain temperatures 
only as a function of the modal radius of particles. Note 
that these are equilibrium temperatures an isolated grain 
would have, if there were no evaporation or grain growth. In 
reality, most of the grains will not survive at temperatures 
higher than the condensation temperature and will quickly 


evaporate (Kimura et al. 2002 ). Also, if gas is present, as 


e.g. in dense brown dwarf atmospheres, grains will exchange 
heat with the gas and both constituents settle to the same 
temperature ( Woitke Helling||2003 ). At much lower den¬ 
sities and in the presence of additional heating sources, like 
e.g. in protoplanetary disc atmospheres, condensates may 
have considerably cooler temperature than the ambient gas 
( [Glassgold, Najita Igea||2004 ). 

Our grain temperatures cover a broad range of distances 
(solid angles). They decline with solid angle. This behaviour 
is not very interesting and, therefore to reveal other more 
interesting features, it is convenient to plot the grain tem¬ 
peratures relative to the temperature of the grey particle, 
j^grey cross-sectiou of grey particle does not depend on 


the frequency. In such case equation (26) gives 


jigvey _ ^ \l /4 


(30) 


This quantity will be used below to normalize and visualize 
the grain temperatures. 


4.5 Radiative accelerations 

Our tables can also be used to calculate the radiative ac¬ 
celeration Ur on a particle of certain modal radius and an 
average mass M®. The radial component of the acceleration 




mil 


CprRdujdu 


(31) 


where c is the speed of light and Cpr is an averaged cross- 
section for radial component of the acceleration given by 




/>' 


/ 


p{c 


pr — j — I OsAoduj ]-^dr, 


(32) 


Al = cos(ai,Ao = cosoo, and ao are angles of the in¬ 
coming and scattered rays with respect to the radial di¬ 
rection, respectively, a is the phase angle of the scattered 
beam (see Fig. |^. The integration is performed over all 
possible directions of the scattered photons doj , over all 
particle sizes dr^ over all directions of incidence du with 
intensity and over all frequencies. The first term is the 
contribution from absorption, which is angle independent 
since the thermal re-emission of the absorbed radiation is 
isotropic. The second and third therms contain the contri¬ 
bution from the scattered photons. Even for close-in planets, 
like WASP 103b, one can assume a point source approxima¬ 
tion and that cosai l,cosao = cos a. Using equations 
@ ,(p^,and ( p^ the following expression is obtained 

Cpr = Ca + (1 — ^)(7s. (33) 


Using the above formula and equations (20) and (27) one 


gets the following expression for the radiative acceleration 
OR = ^ j^s. + {l-g)Cs]B^ (T *)dv (34) 

for the black-body approximation, or 

(jj 




Msttc 


/ 


[C.A-{l-g)Cs]F.{T^)djy 


(35) 


when more precise fluxes, iU, on the surface of the source 
are known. This can be expressed in terms of opacities using 
equations ( 10 ) and 


ttR = - 

TTC 


I 


+ {\-g)atT]FBT*)du. 


(36) 


It is very convenient to express the radiative accelera¬ 
tion as /3 = UR/aG, i.e. relative to the gravitational acceler¬ 
ation: 

GM. 


ac = 




(37) 


Combining equations (36),(37) and eliminating d using uj = 
nRlId^ one gets: 




R 


GM, 


^ J l4T + (1 - g)4T]FAT*)diy. (38) 


This 13 is now independent on the distance or cc. In this way 
we can calculate and tabulate [3 on the same grid of stellar 
effective temperatures T* and particle sizes. We assumed 
M* = Af 0 , R* = Rq and the user will need to adjust our /3 
to his object /3u, using the following relation: 




R^Mq 

R%My 


( 39 ) 
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The numerical integration of the equations was performed 
using a simple trapezoidal rule, as dicussed next. 


4.6 Numerical methods and precision 


The numerical integration of the equations was performed 
using a simple trapezoidal rule on generally non-equidistant 
data points. To assess the robustness of the calculations we 
carried out a number of tests. First, we increased the num¬ 
ber of integration points in the particle size distribution by a 
factor of two. The opacities changed by less than 0.1%, and 
the asymmetry parameter by 0.0002. Second, we extended 
the small size cut-off radius of the particle size distribution 
by a factor of two. The opacities changed by less than 0.2%, 
and the asymmetry parameter by less than 0.0004. Further, 
we extended the upper limit of the particle size distribution 
by a factor of two to larger particles. This is quite an impor¬ 
tant parameter. It resulted in the maximum change in the 
opacity of about 0.3% and in the asymmetry parameter of 
about 0.0003. The scattering opacity is generally much more 
sensitive to this upper cut-off than the absorption opacity 
since, in the Rayleigh scattering regime, it increases with the 
cube of the particle radius. However, the absorption opacity 
may dominate such ’relative (fractional)’ differences for case 
of transparent dust particles with a very low imaginary part 
of the refractive index and low absolute values of absorption 
opacities. 

We also increased the angular resolution (decreased the 
integration step) in the calculations of phase functions by 
a factor of two. The asymmetry parameter changed by less 
than 0.0005 for most of the particle sizes and wavelengths 
and the largest differences in g-values reached about 0.004. 

When extracting fluxes from the model atmosphere ta¬ 
bles we make a simple test and numerically calculate the 
integral of the whole SED: J Fudu. This usually gives a 
small departure from the theoretical value of crT^^. Then, we 
rescale the fluxes to match precisely the theoretical value. 
Consequently, we use the same integration pattern to calcu¬ 
late the dust temperatures and accelerations. 

Our dust opacities have a short wavelength cut-off at 
0.2 micron. For the purpose of calculating the radiative ac¬ 
celerations and dust temperatures we extrapolate them with 
constant boundary value to 0.001 micron. If we were to cut 
them at 0.2 micron this would influence mainly the hottest 
stars, that radiate more at shorter wavelengths, and smaller 
grains that may have significant opacity increase at those 
wavelengths. We found out that the 0.2 micron cut-off did 
not change the radiative accelerations of 0.01 and 100 mi¬ 
cron grains, irradiated by the Tgff = 7000K star, by more 
than 11% and 1%, respectively, and is negligible for cooler 
models. Similarly, our dust opacities have a long wavelength 
cut-off at about 200-500 micron. We have also checked that 
neglecting the opacity at larger wavelengths has a negligible 
effect on the radiative accelerations, since our objects do not 
have a significant flux at those wavelengths. 

However, this long wavelength opacity cut-off may not 
be negligible in the calculations of the dust temperatures. 
If the dust temperatures are low, the integral on the right 


hand side of equation (26) may be underestimated causing 


artificially high dust temperatures. For a 500 micron cut-off 
all grain temperatures that are less than about 30 K will be 
uncertain by more than 2-3%. For ammonia and iron, which 


have about 200 micron opacity cut-off, such uncertainty is 
already noticeable at dust temperatures which are less than 
about 60K. 


5 OUTPUT - DESCRIPTION OF TABLES 

We present our calculations in the form of four separate files 
for each dust species. 

The first file, named e.g. iron_opac_all, deals with opac¬ 
ities. It contains blocks (separated by a blank line) and 
columns. Each block refers to a particular modal particle 
radius Tc- It has five columns: (1) log(rc) in micron, (2) wave¬ 
length in micron, (3) scattering opacity and (4) absorption 
opacity in err? jg (per gram means per gram of iron conden¬ 
sates), and (5) the mean cosine of the phase function. 

The second file, named e.g. iron_phase_all, deals with 
phase functions. It contains several blocks, each block starts 
with log(rc) in micron and refers to a particular modal ra¬ 
dius. Each block has several sub-blocks, which apply to a 
given frequency and are preceded with its value in Hz, the 
mean cosine value, scattering and absorption opacities in 
cw?Ig (per gram means per gram of particular condensates). 
Each sub-block has two columns: phase angle in degrees and 
phase function normalized to An. 

The third file, named e.g. iron.temp, deals with equi¬ 
librium temperatures. It also contains several blocks. Each 
block refers to a particular temperature of the star and 
modal particle radius. There are four columns: stellar ef¬ 
fective temperature [K], log(rc) in micron, solid angle sub¬ 
tended by the star in steradians, and grain equilibrium tem¬ 
perature [K]. 

The fourth file, named e.g. iron_radacc, contains radia¬ 
tive to gravitational acceleration ratio, /3, (3-rd column, di¬ 
mensionless) as a function of the particle size (log(rc) in 
micron, 2-nd column) and the stellar effective temperature 
(in K, 1-st column). It assumes that the mass and radius of 
the irradiating object are Mq,Rq, respectively. 

Opacities were to zero and the phase functions to one in 
the wavelength regions where the refraction index measure¬ 
ments were not available. The provided ’readme’ file contains 
the description of the tables and references to the original 
refractive index measurements of individual species. There 
is also a computer code DISKAVER.E90 in Eortran90 which 
calculates the disc averaged phase functions. 

The behaviour of the phase function, opacity, mean co¬ 
sine, albedo, grain equilibrium temperature, and [3 is cap¬ 
tured in Eigs. |Al] to |A14| 

The dust phase function is plotted using logarithmic 
scale for 1 micron particles. These phase functions indeed 
exhibit strong departures from isotropic scattering with very 
pronounced forward scattering peak. Eor large particles, this 
feature is proportional to the cross-section and rapidly in¬ 
creases with the square of the particle size. Towards low 
X = InajX values (small particles and/or long wavelengths) 
forward scattering decreases and backward scattering in¬ 
creases, the phase functions approach the dipole phase func¬ 
tion, 3(1 + cos(a)/4, and thus have a similar behaviour to 
Thomson or Rayleigh scattering (with some exceptions such 
as iron, see below). The asymmetry parameter reflects the 
behaviour of the phase functions. 

Absorption and scattering opacities as well as the to- 
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tal opacity (absorption plus scattering) and albedo are also 
plotted using the logarithmic scale. Particles relatively small 
compared to the wavelength, with size parameter x << 1 , 
(see e.g. Chap.5 of Bohren & Huffman (1983)) show Rayleigh 
scattering pattern with the scattering efficiency proportional 
to ~ This results in a wavelength dependence and 
a particle-size dependence. Consequently, the scattering 
cross-section is proportional to r® and the opacities per gram 
follow a dependence on particle size. If particles are large 
compared to the wavelength, i.e. x » 1 , Qe approaches a 
value of about 2 and both absorption and scattering opaci¬ 
ties become noticeable grey. 

Albedos are generally high in the optical region. How¬ 
ever, notice that small particles at longer wavelength are 
in the Rayleigh scattering regime, which is not very effec¬ 
tive. As a result, the albedo in this parameter space drops 
by many orders of magnitude to about 10“In the hgures, 
we deliberately saturate this regime and assign a dark brown 
color to albedos smaller than 10“^. The grain temperatures 
decrease with decreasing solid angle. For this reason they are 
plotted relatively to the grey temperatures, (see 

equation (30)) for a solar type star with T* = 5800K on a 
grid of modal particle sizes and solid angles. As the opacities 
of most large particles become progressively more grey this 
results in more grey albedos and equilibrium temperatures. 

/3 values are plotted as a function of particle size for 13 
star/brown dwarf temperatures (see Sec|^. Generally, radia¬ 
tive accelerations can exceed gravity for hotter stars but the 
ratio decreases with decreasing object temperatures. The /3 
values all show a pronounced peak at a certain particle size. 
Assuming a black-body radiation of the object, this peak 
would be at lower grain radii for hotter stars and would 
move towards larger particles for cooler objects. This is due 
to predominant scattering and its dependence for small 
particles together with strong fluxes at the shorter wave¬ 
lengths for hot objects. In the cooler objects, the fluxes move 
to longer wavelengths and so does the Rayleigh regime for 
larger particles, and this enables larger particles to pick-up 
more momentum. We observe that this trend also holds for 
stars without the black-body assumtion, i.e. assuming syn¬ 
thetic spectra from atmospheric models. However, it stalls 
and may even get reversed in objects cooler than about 
1600K. This is associated with the L-T transition. In T- 
type brown dwarfs the opaque silicate dust clouds settle 
into deeper layers, atmosphere cleares, strong molecular ab¬ 
sorption in the infrared region pushes the flux via shorter 
wavelengths reversing the J-K color index. Generally, such 
’non-black-body’ radiative accelerations of larger particles, 
that tend to have more grey opacities, do not differ much 
from the ’black-body’ accelerations, with the exception of 
water and ammonia. However, such ’non-black-body’ accel¬ 
erations of 0.1-0.3 micron particles around cool L,T type 
objects may often be by a factor of 3-5 higher than their 
’black-body’ counterparts. In the next section we briefly 
summarize the most striking features for different species, 
which may be used to identify that particular species in the 
spectra/photometry of exoplanets. 


6 RESULTS 

Alumina/Corundum: (Figs |Al| and |A 2 | absorption opacities 


of smaller particles show a strong 12 micron peak (Koike 


et al. 1995 Zeidler, Posch & Mutschke 2013) which turns 


into a depression for larger particles (Koike et al. 1995 Zei¬ 
dler, Posch Mutschke|20I3 ). Scattering shows a dip at 8-9 


/xm, a peak at 12 fim for smaller particles, and a relatively 
flat behaviour for larger particles. As a consequence we see 
a strong 12 micron feature in the total opacity of small par¬ 
ticles but an enhanced albedo in the 10-30 micron region for 
larger particles. This is accompanied by a sharp drop in the 
mean cosine of the scattering angle g at about 10 fim for 
larger particles. Overall g is positive and approaches unity 
for larger particles and A < 10 gm. One can see strong 
forward scattering but weak backward scattering features. 
Alumina has moderate absorption in the optical regime and 
its equilibrium temperatures approach that of the grey ab¬ 
sorber. Smaller particles, when irradiated by a solar type 
star with solid angle 10“^ — 10“^, tend to have about 50% 
higher temperatures than the grey absorber. Large crystal¬ 
lized corundum particles are quit reflective and have higher 
albedo the IR -far IR regime than 7 alumina. 

Perovskite: (Fig |A3| ) has a strong absorption feature at 
14.1 micron and weaker features at 21.0, 23.8, 33, and 62.3 
micron ( Posch et al.||200^ . Its albedo is very high in the 
near-IR region and remains so for large particles also at the 
longer wavelengths. Signihcant drops in the albedo are ob¬ 
served near the above mentioned absorption features. Mean 
cosine of the phase function reaches negative values at longer 
wavelengths and for larger particles which indicates signih¬ 
cant back-scattering properties. On the other hand, highly 
positive peak g-values are found at II and 20 micron for 
large particles. 

Forsterite: (Fig|^ absorption opacities of forsterite 
show the typical double peak at 10 and 18 micron, which is 
due to stretching and bending of Si-0 molecular bond in sil¬ 
icates. Scattering opacities generally dominate in the visible 
region and become progressively more important at longer 
wavelengths as the particle size increases. Smaller grain sizes 
exhibit the A“^ dependence typical of the Rayleigh regime. 
These opacity properties are well known and we will now fo¬ 
cus on the behaviour of the phase functions. The mean cosine 
of the scattering angle is positive, which indicates predomi¬ 
nance of the forward scattering. Forward and backward scat¬ 
tering increase towards shorter wavelengths or larger parti¬ 
cles. As the particle size increases backward scattering first 
disappears at shorter wavelength and then reemerges and 
expands to progressively longer wavelength. Phase functions 
also show strong wavelength dependence in the vicinity of 
10-18 gm features. Forsterite absorbs little flux in the opti¬ 
cal region which results in its equilibrium temperatures for 
solar type stars being signihcantly lower than those of the 
grey absorber. Its albedo is very close to unity in the optical 
and NIR regions but large particles remain reflective also in 
the far-IR and sub-mm region. Broad absorption bands at 
10-20 micron signihcantly reduce its albedo. 

Olivine: (Fig |A5| ) with 50% of iron has considerably 
stronger absorption in UV, visible, and NIR region than 
forsterite for smaller and medium size particles. Small 
olivine particles have slightly higher total opacity in the 
NIR region. The 10 micron feature is sharper, especially its 
blue edge, while the 20 micron feature is slightly broader. 
The albedo of olivine is significantly smaller than that of 
forsterite in the whole optical and NIR region. Gonsequently, 
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the equilibrium grain temperatures of medium to small par¬ 
ticles are several times higher than those of forsterite. Larger 
olivine particles have temperatures similar to the grey par¬ 
ticles. This is a consequence of the imaginary part of the re¬ 
fractive index of olivine, which is much higher in the optical 
and NIR regions. Radiative accelerations of the small olivine 
particles are significantly higher than those of forsterite and 
the maximum is less pronounced. The forward scattering of 1 
micron particles in the optical region is narrower while back 
scattering peak diminished and, consequently, the asymme¬ 
try parameter of medium size grains is larger. 

Enstatite: (Fig@ behaves in a similar way as 
forsterite. Its absorption opacity in the visible and near in¬ 
frared regions is small, even smaller than that of forsterite. 
Consequently, its equilibrium temperatures are consider¬ 
ably below the grey absorber and even lower than those 
of forsterite. The phase functions and mean cosine of the 
phase angle show even more pronounced wavelength depen¬ 
dence in the 4-20 iim region. The albedo of enstatite has 
similar behaviour to that of forsterite. Both iron-free sili¬ 
cates, forsterite and enstatite, show a pronounced peak in /3 
values with relatively low radiative accelerations, especially 
for small grains, when compared to other species. 

Pyroxenes: (Figs |A7] and |A8| ) even a small amount of 
iron considerably raises the imaginary part of the refractive 


index in the UV, visible, and near IR regions (Dorschner 
et al.||1995 ). This raises the absorption opacity and equilib¬ 


rium grain temperatures. Particularly small grains, at rela¬ 
tively close distance to a solar type star, with a 20% amount 
of iron, may become almost 2 times hotter relative to the 
grey absorber. Larger particles can still remain relatively 
cool at larger distances. Its albedo shows similarities to iron- 
free silicates but we start to see an enhanced absorption in 
the UV for larger particles and accompanied reduction of the 
albedo. This trend gets even more pronounced and spreads 
to longer wavelengths and smaller particles if the amount of 
iron is increased to 60%. The presence of iron has also an 
effect on the phase functions and increases the mean cosine 
of the scattering angle for A < 4 fim. Backward scattering is 
apparent in visible and near IR regions, however, it dimin¬ 
ishes for large particles in UV region and this trend is more 
pronounced with higher iron content. 

Iron: (Fig. |A9[ ) contrary to silicates, iron has consider¬ 
able absorption opacity in the visible and IR, which gen¬ 
erally smoothly declines towards longer wavelengths and 
larger particles. Its absorption and scattering opacities are 
almost featureless (grey) apart from the Rayleigh scattering 
regime for smaller particles. The mean cosine of the scat¬ 
tering angle increases with the particle size and towards 
shorter wavelenths. It may have significant negative val¬ 
ues, especially at longer wavelengths, which is due to a very 
pronounced backward scattering peak. For larger particles, 
forward scattering overcomes the backward scattering at 
shorter and progressively at longer wavelengths. These fea¬ 
tures can be used to identify iron in exoplanets. Highly effec¬ 
tive absorption properties in the optical and IR regions heat 
iron grains to significantly higher temperatures as expected 
for the grey absorber. Especially, small iron grains at large 
distances from solar-type stars can get heated to tempera¬ 
tures more than a factor of 3 higher than the grey absorber. 
As a result of its opacity behaviour, large iron particles are 
efficient at scattering in the IR and far IR regime. However, 


compared to silicates, iron particles (especially small ones) 
have quite low albedo in the UV and -visible region. /3 values 
peak at slightly smaller grains than most of other species. 
Small iron grains will acquire significantly higher radiative 
accelerations than other species when irradiated by stars. 
This is due to high iron opacities in the optical as well as 
strong back-scattering properties which push the particles 
away from the source. On the contrary, larger iron grains 
tend to have slightly smaller /3 values for all source temper¬ 
atures. 


Carbon: (Figs. AlO and All) differs from most of the 
species but some features resemble those of pure iron. Car¬ 
bon at higher temperatures and higher graphitization, has a 
total opacity which smoothly declines towards longer wave¬ 
lengths and larger particles. Apart from the slope, the dom¬ 
inant ’feature’ is the fact that it is almost featureless. There 
is a drop in the absorption and scattering opacities of small 
particles in the UV. The most striking characteristic is that 
its albedo in the optical and IR regions is much smaller 
than that of other species, even smaller than that of iron. 
However, the albedo of large particles at far-IR and sub-mm 
wavelengths is approaching unity. The asymmetry parame¬ 
ter is mostly positive and smoothly increases towards the 
shorter wavelength. However, it becomes negative for larger 
particles in the sub-mm region, a feature indicating stronger 
backward scattering but not as pronounced as in iron. The 
phase functions lack any abrupt features too. 

Carbon at cooler temperatures has a steeper gradient 
of the total opacity in the IR region but more shallow in 


the far-IR and sub-mm region (Jager, Mutschke V Henning 


1998). Its albedo features an enhanced band traveling from 


medium size particles at NIR region towards larger particles 
at sub-mm region caused by an enhanced scattering at the 
blue edge of the Rayleigh scattering regime. Apart from this 
feature, the opacities are almost perfectly grey and feature¬ 
less. Contrary to its behaviour at hotter temperatures, the 
asymmetry parameter is everywhere positive and increasing 
smoothly towards larger particles and shorter wavelengths. 
Dust temperatures of larger particles and closer to the solar 
type star are slightly cooler while smaller particles which are 
farther from the star are significantly hotter than the grey 
absorber. Small carbon grains acquire the highest radiative 
accelerations of all the species studied here. 

Water iee: (Fig. |A12 ) absorbs very little in the visible 
but has strong absorption and scattering features in the IR, 
particularly at about 3 micron, and shallower features at 
about 12 and 40-50 jam. This is accompanied by the vari¬ 
ability of phase functions at these wavelengths. Mean cosine 
values of larger particles show extra particle size sensitiv¬ 
ity in the 3-10 lam region and is larger than zero. However, 
broad back scattering feature is present in the visible region. 
Its extremely low absorption in the optical causes the grains 
to be significantly cooler than those of the grey absorber for 
hotter stars. Its high absorption in the IR, however, might 
heat smaller grains at large distances to become hotter than 
grey absorber when irradiated by cooler brown dwarfs. The 
albedo of water ice is very high, essentially 1, for all particles 
from UV to NIR. Large ice particles are slightly more reflec¬ 
tive than liquid drops in the far IR and sub-mm regimes. In 
the IR one can see the water absorption bands which lower 
its albedo. 

Water liquid: (Fig. |A13 ) has similar features as ice. It 
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has slightly more pronounced absorption features, e.g. a nib 
at about 6 /xm, and slightly higher UV and blue absorption. 
Mean cosine values of larger particles show enhanced particle 
size sensitivity in the 3-10 fim region. 

Ammonia: (Fig. A14) has very sharp and strong absorp¬ 
tion and scattering features at about 3, 9.4, 26, 70, 110-140 
jam, which may interfere with the watter and silicate fea¬ 
tures. Its absorption in the visible is much stronger than that 
of watter. Apart from forward scattering, its phase functions 
exhibit also back scattering at the shorter wavelengths and 
strong variability with wavelength at the above mentioned 
spectral regions which can help to identify ammonia in re¬ 
flected light. The equilibrium grain temperatures of ammo¬ 
nia are similar to those of water but larger particles tend to 
follow the trend of the grey absorber while smaller particles 
at medium distances are even cooler. Overall, the albedo of 
ammonia has similar characteristics as that of water ice. A 
number of absorption features can be seen which reduce its 
albedo in the IR and far IR regions. Medium to large ammo¬ 
nia grains also experience strong radiative push compared 
to other particles. 


7 CONCLUSIONS 

We have calculated phase functions, asymmetry parameters, 
absorption and scattering opacities, single scattering albe¬ 
dos, equilibrium dust grain temperatures, and radiative ac¬ 
celerations for about a dozen of the species most commonly 
encountered in the extrasolar planet environment. Mie the¬ 
ory and spherical homogeneous grains are assumed. The op¬ 
tical properties are integrated across a poly-dispersed Deir- 
mendjian particle size distribution with different modal par¬ 
ticle radii. These distributions are relatively narrow. Conse¬ 
quently, they might be understood as single sized particles 
with the radius equal to the modal radius or used to as¬ 
semble and study an optional broader particle size distri¬ 
bution, if necessary. These dust properties are presented in 
the form of tables as a function of wavelength and modal 
particle size. Equilibrium grain temperatures are a function 
of modal particle size, host star effective temperature, and 
solid angle subtended by the star. Radiative accelerations 
are tabulated as a function of modal particle radius and 
host star temperatures. Dust opacities are deliberately cal¬ 
culated per gram of the dust material so that the user has 
more flexibility and can place these dust particles into any 
environment independently on the gas to dust ratio. Tables 
will be maintained, potentially expanded, and available for 
download from: 

\protect\vrule widthOpt\protect\href{http://www.ta3. 

There is also a code which allows the user to take into ac¬ 
count, an optional, finite dimension of the source of light 
(star) and calculate grid of disc-averaged phase functions. 
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Figure Al. Alumina. Very-top: absorption and scattering opacities. Top left: total opacity (scattering plus absorption). Top right: single 
scattering albedo. Middle left: mean cosine of the scattering angle. Middle right: phase function for modal particle radius of 1 micron. 
Bottom left: Grain temperature relative to the grey absorber for a solar type star with T* = SSOOAT. Bottom right: Radiative acceleration 
relative to the gravity of the central object for 13 stellar temperatures, from the top: 7000, 5800, 5000, 4500, 4000, 3500, 3000, 2500, 
2000, 1600, 1200, 900 and 700 K. 





















































log par. size [mic] log par. size [mic] log par. size [mic] 


Tables of dust properties in exoplanets 15 


log absorption opacity [cm^2/g] 



0.2 1 10 100 500 

lannbda [mic] 


log opacity [cm^2/g] 



0.2 1 10 100 500 

lambda [mic] 



1 

0.8 

0.6 

0.4 

0.2 

0 

- 0.2 

-0.4 


0.2 1 


10 100 


500 



0.2 1 10 100 500 

lambda [mic] 



0.2 1 10 100 500 

lambda [mic] 



lambda [mic] 


lambda [mic] 


Figure A2. Corundum. Pictures and notation are the same as in the previous figure. Dust temperatures and accelerations are not 
calculated in this case because the opacities do not cover the whole spectral region. 
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Figure A3. Perovskite. Pictures and notation are the same as in the previous figure. Dust temperatures and accelerations are not 
calculated in this case because the opacities do not cover the whole spectral region. 
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Figure A4. Forsterite. Pictures and notation are the same as in the previous figure. 
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Figure A5. Olivine with 50% of iron. Pictures and notation are the same as in the previous figure. 
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Figure A6. Enstatite. Pictures and notation are the same as in the previous figure. 
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Figure A7. Pyroxene with 20% or iron. Pictures and notation are the same as in the previous figure. 
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Figure A8. Pyroxene with 60% or iron. Pictures and notation are the same as in the previous figure. 
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Figure A9. Iron. Pictures and notation are the same as in the previous figure. 
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Figure AlO. Carbon assuming optical data for 1000 C. Pictures and notation are the same as in the previous figure. 
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Figure All. Carbon assuming optical data for 400 C. Pictures and notation are the same as in the previous figure. 
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Figure A12. Water ice. Pictures and notation are the same as in the previous figure. 
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Figure A13. Liquid water. Pictures and notation are the same as in the previous figure. 
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Figure A14. Ammonia. Pictures and notation are the same as in the previous figure. 



























































